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Abstract. We obtain a sharp local well-posedness result for the Gradient 
Nonlinear Wave Equation on a nonsmooth curved background. In the process 
we introduce variable coefficient versions of Bourgain's X s,b spaces, and use a 
trilinear multiscale wave packet decomposition in order to prove a key trilinear 
estimate. 



1. Introduction 

In this article we are investigating the issue of local well-posedness for a variable 
coefficient semilinear wave equation in 4 + 1 dimensions. To describe the context 
and motivate the interest in our problem we introduce three related equations. We 
begin with a generic gradient NLW equation in R ,l+1 , 

(1) Du = T(u)(Vuf 
with the nonlincarity 

r(«)(Vu) 2 = q ij (u)diud jU 
where q lJ are smooth functions and the standard summation convention is used. 
Then we move on to a similar equation but on a curved background, 

(2) D g u = T(u)(Vw) 2 

with n\ g — g VJ didj, where the summation occurs from to n and the index 
stands for the time variable. To insure hyperbolicity we assume that the matrix g lJ 
has signature (l,n) and the time level sets xq = const are space-like, i.e. g oa > 0. 
In effect to simplify some of the computations we make the harmless assumption 

Finally, we consider a corresponding quasilinear equation 

(3) n g{u) u = T(u)(Vu) 2 

with similar assumptions on the matrix g. 

In all three cases we are interested in the local well-posedness of the Cauchy 
problem in Sobolev spaces H s (R n ) x H s ~ 1 (W l ) with initial data 

(4) u(0,x) — uq(x), dtu(0,x) = ui(x) 

The first equation (fT]) is the best understood so far, and is known to be locally 
well-posed for s in the range 

c n n + 5 
s > max{-, — -p } 

This range is sharp. The ^ obstruction comes from scaling, while the is related 
to concentration along light rays, see Lindblad [9 . The proof of the positive result 
is fairly straightforward in dimension 2 + 1 and 3 + 1, where it suffices to rely on the 
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Strichartz estimates. In 4 + 1 dimensions this no longer works and one needs to use 
instead the X s,e spaces, see Foschi-Klainerman [3]. These are multiplier weighted 
L 2 spaces associated to the wave operator as the Sobolev spaces H s are connected 
to the Laplace operator A, see Klainerman-Machedon [5]: 

(5) IMI*... = ||(1 + iCl 2 ) 4 ' (1 + ||r| - |£|| 2 )^ • |fi(r,OIII*. 

where u = u(r, £) is the space-time Fourier transform of function u = u(t,x). 
Finally, in the most difficult case, n > 5, this was proved by Tataru |14j . using a 
suitable modification of the X s ' spaces, needed in order to control the interaction 
of high and low frequencies in the multiplicative estimates. 

For the quasilinear problem ([3]) the sharp result is only known to hold in di- 
mensions n — 2, 3. This was proved by Smith- Tataru [12] (see also Lindblad's 
counterexample |10jV The argument there still requires the use of Strichartz esti- 
mates. These are derived from a wave packet parametrix construction for a wave 
equation with very rough coefficients, which in turn is obtained via a very deli- 
cate analysis of the Hamilton flow. A different proof of this result in the special 
case of the Einstein vacuum equation was independently obtained by Klainerman- 
Rodnianski [6], [8], [7]. In dimensions n > 4 it is still unclear which is the optimal 
threshold, the best results so far being contained in the above mentioned paper of 
Smith- Tataru |12j and in an earlier one, Tataru |15j : 

n = 4, 5 

n > 6 

In the same direction but somewhat closer in spirit to the present paper is Bahouri 
and Chcmin's work [2][I]. The equation considered there is still quasilinear, but 
the main estimates are frequency localized versions of the Strichartz estimates for 
the wave equation on a rough background. 

As an intermediate step toward understanding the higher dimensional quasilinear 
problem, we consider here the semilinear problem on a curved background and we 
prove the sharp result: 

Theorem 1.1. Let n = 4 and assume that the coefficients g 1 ^ satisfy d 2 g £ L 2 L°° . 
Then the Cauchy problem @, ([U is locally well-posed in H s x H"~ 1 for s > |. 

Here well-posedness is understood in the strongest sense, i.e. the solutions have 
Lipschitz dependence on the initial data and they exist on a time interval which 
only depends on the size of the initial data. 

One contribution of the present paper is to introduce variable coefficient versions 
of the X s ' b spaces, study their properties and obtain the corresponding Strichartz 
type embeddings. However, the main novelty, contained in the last two sections, 
is a new method, based on a trilinear wave packet decomposition, to prove a key 
trilinear bound which cannot be obtained directly from the Strichartz estimates. 

The first step in the proof is to reduce the problem to the case when the initial 
data is small, using scaling and the finite speed of propagation. This is a routine 
argument for which we refer the reader to [12] . Once we know that the initial data 
is small, we can fix the time interval and set it to [—1,1]. This will be the case 
throughout the rest of the paper. 
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To solve the problem for small data we use a fixed point argument. Let S(uo, Ui) 
and D" 1 be respectively the homogeneous and inhomogeneous solution operators 

(6) n g S(uo,ui) = 0, S(u o ,ui)(0) = uq, dtS(uo,ui)(0) = u\ 

(7) a^H) = h, (nfHM = o, dtinfHM = o 

Then a solution u for $2$ in [—1,1] is also a fixed point for the functional 

(8) F(u) = S(u , ui) + □; 1 (I»(V U ) 2 ) 

In order to apply a fixed point argument for F we need to find two Banach spaces 
X and Y for which the following mapping properties hold: 

(9) ||5(u ,Ul)||x < ||(W0,Ml)||iTa X H»-l 

(10) lin-^Hx < \\H\\ Y 

(11) ||u • w\\x < IM|x|MU 

(12) l|r(n)||x<C(|H| ic »)(l + ||u|| 5 x ) 

(13) \\u ■ w\\ Y < IMUHHIy 

(14) ||Vu- Viy||y < \\v\\ x ■ \\w\\x 

where C = COIuHlcxj) is a constant that depends solely on ||w||i,oo. In the flat case 
([T]), for dimension n = 4, one can make this argument work by choosing 

X = X s ' 9 Y = X s - 1 ' - 1 

with 

(15) a = 9 + 1 9 > \ 

For our problem the challenge is twofold: first we need to find suitable variable 
coefficient versions for the X s ' spaces and then, in this new context, prove the 
corresponding estimates (|9|- ((l"4"|) . 

Such spaces were previously introduced by Tataru |13j . where they are used in 
the context of a unique continuation problem. There, for a hyperbolic operator P 
one defines 

X s,o = H s^ X s,i = | u g h s \Pu e H 3 - 1 } 

Then all the other spaces are defined through interpolation and duality. 

In this article we choose to follow a different path based on dyadic decompositions 
with respect to the spatial frequency and the distance to the characteristic cone. 
Likely one should be able to prove that the two approaches are equivalent, but we 
choose not to pursue this here. 

Our article is structured as follows. In the next section we define the X s ' 
spaces and prove that they satisfy the linear estimates ([9|), (flUf . Our definition of 
the X s ' 6 is slightly different from the standard one (|5|) in the constant coefficient 
case. Precisely, in the constant coefficient case our definition gives 

(16) ||«|| x ... w ||(l + |e| 2 )*(l+||T|-|e|| 2 )^w(T,e)|| L2 + ||D U || i? ^ +e - 2 , < 9 < 1 

and one can see that the second term above alters the behavior at high modulations 
|t| 3> |£|. Correspondingly, for negative 9 we have 

(17) ||u|| xs , e «||(l + |e| 2 )5(l + IM-|e|| 2 )^u(r,0llL 2 + hll L? ^, -l<0<0 
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This change is consistent with scaling and simplifies somewhat the study of high 
modulation interactions. 

In Section [3] we discuss the Strichartz estimates for O g , which translate into 
embeddings for the X s ' spaces. These turn out to suffice for the proof of the 
algebra properties (fTTj) - (fT5|) and for the high-high frequency interactions in p4[) . 

The difficult part is to study the high- low frequency interactions in ([33]). For 
this we first take advantage of the duality relation 

(18) (X s ' 9 + L 2 H s+e Y = X- s >- e s€R, 0<6»<i 

This is consistent with (fT6|) and (fT"7)) . Using this duality, after factoring out high 
modulation interactions, the bound (|14[) is transformed into the trilinear estimate: 



(19) 



u • v ■ w dxdt 



< ||u|| x i-.,i-e||w|| x .-i,e|H| x . 



with (s,0) verifying (fT5|) . The last section of the paper is devoted to proving this 
bound. The argument is based on a multiscale trilinear wave packet decomposition 
for linear waves. 

2. The X s ' 6 spaces 

We first introduce Littlewood-Paley decompositions. As a general rule, all fre- 
quency localizations in the sequel are only with respect to the spatial variables. 
There is a single exception to this. Precisely, the coefficients g* 3 are truncated 
using space-time multipliers. In order for these truncations to work, we need for 
these coefficients to be defined globally in time. Hence we assume they have been 
extended to functions with similar properties in all of R n+1 . 

Let 4> be a smooth function supported in {| < |£| < 2} with the property that 

oo 

1= # 2-j o 

j=—oo 

We consider a spatial Littlewood-Paley decomposition, 

oo 



where for dyadic A > 1 we have 
while Si incorporates the low frequency contribution in {|£| < 1}. Set 



Sx (0 



We will also use spatial multipliers S\ with slightly larger support, with S\S\ = S\. 
We say that a function u is localized at frequency A if its Fourier transform is 
supported in the annulus {j < |£| < 8A}. 

For the paradifferential type calculus we also need to truncate the coefficients of 
\3 g in frequency. Given D g in ^fy we define the modified operators 
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In the sequel we omit the space and time variables in our function space notations, 
i.e. LP := L p x t , L 2 H S := L 2 H S X , LPL q := L p t L%, etc. We are ready now to define 
our spaces: 

Definition 2.1. Let 6 6 (0,1) and s € M. Then X s ' 8 is the space of functions 
u € L 2 (— 1, 1; H S (W 1 )) for which the following norm is finite: 



(20) 



{oo A oo A ^ 

EEii UA > d iix^ ; "^EE^^ 
A=l d=l A=l d=l ) 

where X, d take dyadic values and 

(2i) IIumII' ... = x 2s d 26 \\u x , d \\h + \ 2s - 2 d 26 - 2 \\n g n x , d \\l, 

We also define the space of functions for which the following norm is 

finite: 



(22) 1 



*A 

A=l d=l 



oo A *j 
A=ld=l J 



Remark 2.2. Intuitively d stands for the modulation of the u\.d piece. Indeed, in 
the constant coefficient case one can easily see that u x ^d mainly contributes to u 
in the region where \\t\ — |£|| w d. The condition 1 < d is related to the spatial 
localization on the unit scale in our problem. The condition d < X reflects the fact 
that at high modulation we use a simpler structure, see e.g. (|16p . (|17[) . 

Remark 2.3. The cutoff at frequency less than y/X for the coefficients d g is related 
to the regularity of the coefficients, d 2 g € L 2 L°° . This implies that O g>t/J .u Xt d is 
an allowable error term. 



We begin our analysis of the X s * spaces with a simple observation, namely that 
without any restriction in generality one can assume that the functions and 
fx,d in Definition 12. II are localized at frequency A. Precisely, we have the stronger 
result: 

Lemma 2.4. The following estimate holds: 

(23) A*-V||VSa«|U« + \ a - 1 ep- 1 \\n g<sn .Sxv\\ V i < \\v\\ x .,» 

Proof. We first bound the time derivatives of v in negative Sobolev spaces, 

(24) A s d e (||a t 2 «|| i 2 (H - 2+A 2 L2) + \\d t v\\ L 2 {H -i +XL2) ) < \\v\\ x s,o 

A ,d 

This follows by Cauchy-Schwartz from the interpolation inequality 

ll<9t u lll,= (ff-i+A£2) < {\\d 2 v\\ L 2 (H -2 +X 2 L 2 ) + 

combined with the bound 

\\8$v\\ L n H -, +x »v>) < A- a ||D fl v|| £ a + \\d t v\\ L 2 {H -i +XL 2 } + \\v\\ L 2 
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To prove this last estimate we only use the L°° regularity of g together with the 
condition g 00 = 1. Then we need the fixed time bounds 

\\9 < ^xd x d t v\\ H -2 +X 2 L 2 < \\d t v\\ H -i +XL 2 

\\9 < -/xdlv\\ H -2 +X 2 L 2 < \\v\\ L 2 
They are similar, so we only discuss the second one. We write 

and use the uniform bounds 

\9<Vx\ S 1> l^<Vll £ A > I^<VaI £ A ' 
This concludes the proof of (fM|) . 

The first term in (f25| is directly bounded using For the second it suffices 

to prove the commutator estimate 

(25) \\[D g< ^S x ]v\\ L 2<X\\v\\ L 2 + \\d t v\\ L 2 

We have 

[□ 9<vx , s A ] = [ 5<vx , s a ]sa + [ 5<vx , s A ]^ 

and the commutators are localized at frequency A so the spatial derivatives only 
contribute factors of A. Hence (|25[) follows from the standard commutator estimate 

ll[<?<v^ A ]|U^ L2 <\- X \\V9\\l~ 

□ 

Applying the above Lemma with S x replaced by S x we obtain 
Corollary 2.5. One can replace the X^' e d norm in the definition of X s ' and 

jp-l,0-l fry n0rm 



= A a - 1 d°||V«||La + A'- 1 d fl - 1 ||n 



For the proof of the duality relation (fT5|) it is convenient to work with a selfadjoint 
operator. Thus we consider the selfadjoint counterpart O g of O g 

n g = ih.r'i) J 

Then for v localized at frequency A we commute and estimate the frequency local- 
ized difference 

\p g< ^v-a g< ^v\\ L 2<\\wv\\ L 2 

This leads directly to 



Corollary 2.6. One can replace the D g< _ operator in the definition of X s ' 9 and 
X s-1, 1 by the similar operator in divergence form S< ^_. 

As a consequence of the second part of (p?5|) we have 

Corollary 2.7. The following embedding holds for —1 < 8 < 0: 

X s,9 c L 2 H s+9 
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Another use of this is to establish energy estimates. A direct application of 
energy estimates for the wave equation yields the bound 

l|V«|||„ i2 <||V«|| 2 L2 + ||V«|| i2 ||D s i;|| L2 

This leads to 

(26) \'- 1 t?-i\\VSxv\\i,- L 2<\\v\\ jt ..e 

A, d 

Going back to Definition 12. H this implies 
Corollary 2.8. Assume that 8 > |. Then 

(27) ||u||i-ir. + KHi-jr-i £ ||«IU».« 



To prove the estimates © and (JTUJ) in the context of the X s,e spaces we need 
to switch from the frequency truncated coefficients to the full coefficients <? y . The 
tool needed to do that is contained in the following: 

Lemma 2.9. Assume that < s < 3. Then the following fixed time estimate holds: 

oo 

(28) EA^H&Cffwx* ~ (^(II^IU»)) 2 ||«ll^- 2 

A=l 

where M stands for the maximal function with respect to time. We also have the 
dual estimate 

oo oo 

(29) llEfwx&Allfl*-. <£ A2(1 ~ s) !IMIi 2 

A=l A=l 

Proof. We take a Littlewood-Paley decomposition of both factors, 

oo oo 

The (/i, i/) term is nonzero only in the following situations: 

(i) v <C A, /i ~ A. Then we estimate 

II&Gm^iu. < IMMKlUa < v- 2 M{\\d 2 g \\ L ~)\\uj\ L * 

and use the square summability with respect to A together with the relation v~ 2 < 
A" 1 . 

(ii) v ss A, fi <C A. Then 

WSxiSvuJU* < IMU-KIU 2 < \- 2 M(\\d 2 g\\ L -)\\u^\\ L i 

This is tight only when s = 3 and fi = 1, otherwise there is a gain which insures 
the summability in A, fi. 

(iii) v fts fj, > A. Then 

I|Sa(<M»„)||j? £ \\9v\\l~\K\\ L 2 < ^ 2 M(\\d 2 g\\ L ~)\\u^\\ L2 
This is always stronger than we need. The proof of the lemma is concluded. □ 
We now establish some simple properties of the linear equation 

(30) D g u = f, u(0)=tto, ut(0)=ui. 
Then 
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Lemma 2.10. The linear equation ([30]) is well-posed in H s x H s 1 for < s < 3. 

The proof follows easily from energy estimates, see [16]. 
We use this to prove |9j), namely 

Lemma 2.11. Assume that < s < 3 and 9 > 0. T/ien the solution u to (|30p 
weri/ies 

IMIx- J$ ||«0||ff. + IMU-l + ll/IU>H-i 

Proof. We decompose the solution it as 

oo 

it = S\S\u 

A=l 



and think of this as a part of the sum in (|20|) which corresponds to d = 1 . Then 

Ni^..<f;ii5 A «ii^,. 



A=l 



« j> 2 i&u|li. + \ 2 ^\\n g<v - x s xu \\l 2 

oo 

< h\\l* H . + J2 x2{s ~ 1) \\ a 9< v-Au-s x a g u\\ 2 L2 + ll/ll 2 ^-, 

A=l 

The first term is easily controlled by energy estimates. The second is decomposed 
as follows: 

For the commutator we use the fixed time bound (|25|) along with square summa- 
bility in A. The second part is controlled by (|28|) . 

□ 

The result in the next Lemma implies the estimate (|TU| for the spaces X, Y: 

Lemma 2.12. Assume that < s < 3 and | < 6 < 1. Then the operator D^ 1 has 
the mapping property 

n- 1 : X s - 1 ' 9 " 1 -+ X s ' 6 



Proof. Let / £ X s ■ . We use the representation in 

oo A 
A=l d=l 

By Definition (|2.ip the function 

oo A 
A=l d=l 

belongs to X s ' 9 . The difference v = u — O^ 1 / solves 

n gV = n g u-f, «(o) = t»(o), «t(o) = ut(o) 
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To estimate it we use Lemma 12.111 The initial data is controlled due to Corol- 
larv !2.8l so it remains to bound the inhomogeneous term in L?H S ~ X . Thus we need 
to show that 

oo A 

(31) || E n^SxhAl'H-i Z E II/mII^. s 

X=ld=l X,d 

Considering the trace regularity result in Corollary 12.81 this would follow from 

oo A 

IIE^vx^/aIIW-^EH^II!^.-. A = £/m 

A=l A d=l 

which in turn is a consequence of the fixed time bound (|29p . 

□ 

We finish this section by proving a key duality relation between X s,e spaces with 
positive, respectively negative 9. 

Lemma 2.13. For < 9 < = we have the duality relation 

(32) X- s '-° = {X sfi + L 2 H s+e )' 
Proof, a) We first show that 



C (X s ' 8 + L z H s+a ) 



2 tjs+8\I 



From Corollary 1 2. 71 we obtain X s - 6 C (L 2 H s+e )' . It remains to prove the bound 



u ■ f dxdt 



< 



WAx^o WfWx— .- 



We consider Littlewood-Paley decompositions of u and v as in Definition 12. 1[ 

oo A oo A 

m = EE Sau «' / = /o + EE s f<vA s '^^.rf 



A=l d=l 



X=l d=l 



with D g ^_ in divergence form, see Corollary 12.61 The summation with respect 
to A is essentially diagonal therefore it follows by orthogonality. To handle the d 
summation it suffices to obtain the off-diagonal decay 



Sxux,d ± ■ U g<v -S\f\.d 2 dxdt 



<i 



a-; 



d-2 

I/a.^II 



If c?2 < d\ then this follows directly from (f2"Tj) and (j!?3"|) . Otherwise we integrate by 
parts 



S\u\,d 1 -O g<v -.S\f\,d 2 dxdt= / \3 g< ^S\u\ l d 1 • S\fx,d 2 dxdt 



+ J (Sxux,di • 9 < ^/xd a S\f\ t d 2 — g < ^/ x -daSxux,di • S\f\,d 2 )dx 

For the first term we use (|23|) and (|2Tj) . For the second we use the trace regularity 
result in (|26| . 

b) We now show that 

{X s > e +L 2 H s+e )' C x- s ~ e 
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Let T be a bounded linear functional on X s ' 6 + L 2 H s+e . Due to the second term 
we can identify T with a function u € L 2 H~ s ~ e . 

On the other hand, we can apply it to functions v S X s,e of the form 



oo A 



y = ee sxvxm 



\=i d=i 

Then we must have the bound 



Tv\* < \\v\\\„ e < £ IK4^ S E (>? s d 2e \\v.A\h + A 2s - 2 rf 2e - 2 ||n fl<vT «A )rf ||i 3 ) 



A,d X.d 



Given the definition of the X^ d norms, using succesively the Hahn-Banach theorem 
and Riesz's theorem it follows that we can find functions fx,d an d h\d with 

oo A 

( 33 ) EE^^H/^II^ + A^-^^-^H^IIi. = M < oo 

A=l d=l 

so that 

oo A p 

Tv = ^2 E / f X d Vx > d + hx > d ' [ - 1 9<vt Wa ^ dxdt 

A=l d=l 

In particular this must hold for v of the form v — S\V\ d, 



U S\v x .ddxdt = J fx.d v\, d + h\,d ■ D g< ^v x , d dxdt 
For each A, d this yields 

S\u = /A,d + U g< ^-h\ id 
Then we can represent S\u in the form 

A 

2 

(34) S\u = f x ,i + ^2 n g< ^-u\, d + □ 5<v _/ia,a Ux,d = hx,d-hx 

d=l 

This yields for u the representation 

OO / 2 

(35) « = E ^ KA,i + E fi s<vx u M + fi ff <^^,A 

A=l \ d=l 

This is very close to but not exactly the form in (|2"2")l . However the multipliers 
S\ can be easily replaced by 5a by reapplying the Paley-Littlewood decomposition 
on the right, and then S\ can be commuted to the right of n s ^_ due to the 
Corollary 12 . 51 and the commutator bound (|25[) . Hence we have 

A/2 



M\ 2 x -„-e < E A_2S |l/A,l|||2 + E ll«A,«l|lii-.,i-. + ||/ia,a||^x-»,x- 
A=l \ d=l 



and due to (|33[) it remains to bound the right hand side by 

M+\\u\\ 2 L2H ^e 
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There is nothing to do for the /a,i term. On the other hand we can bound 

IK*|&i-..i-. < A 2 ^ 1_s ^d 2 ^ 1_ ^||ttA,<i||i2 + \~ 2s d~ 2e \\n g<vx u\ t d\\%2 

= \W->W l - e )\\h Xld - hx&wl* + x~ 2s d- 2e \\h, d - f x , 2d \\h 

<\- 2S d- 2e (\\h, d \\h + \\fx,2 d \\h) 

+ A 2(l- s ) d 2(l- e)(| |^ 2d| |2 2 + || ^ d| |2 2) 

Finally, for the last term we have 

HfrA.All^-..!-. < x 2{1 - s) x 2{1 - 9) \\h XA \\ 2 L2 + x- 2s x- 2S \\d g ,h x , x \\ 2 L2 

A , A ^ v 

= \ 2il - s) \ 2il - e) \\h x , x \\h + X- 2s X- 2e \\S x u - f x , x \\ 2 L2 

< A 2(1- S ) A 2(1-0) ||/IA)A|| 2 2 + A -2. A -M|| /A a ||2 2 + || 5aU || ^ 2ff _ s _ e 

The proof is concluded. □ 

3. STRICHARTZ ESTIMATES AND APPLICATIONS. 

The Strichartz estimates for the variable coefficient wave equation, as proved in 
[15j . have the form: 

Theorem 3.1. (Tataru |15j ) Assume that the coefficients g 13 ofO g satisfy d 2 g lJ g 
i 1 ^ 00 . Then the solutions to the wave equation in n + 1 dimensions satisfy the 
bounds 

(36) H^VuIIlpl, < ||u(0)|| H i + ||ttt(0)|| £a + \\a g u\\ L i L 2 

where 

n 1 n 2 n — 1 n — 1 

37 o- =-- + - + -, - + < — — , 2<p<oo, 2<q<oo 

2 P 9 P 9 2 

Applying this bound on an interval I of size e 2 we obtain by Cauchy-Schwartz 

\\D a Vu\\ LV(I . Lq) < ~\\u\\ H i( IxMv .j + e\\a g u\\ L 2 {IxRn) , e < 1 

Summing up over small intervals this extends to intervals of arbitrary lengths. 
Optimizing over e yields 

(38) \\D°Vu\\ 2 LPLq < \\u\\ 2 H1 + || u || ff i|p 9 u|| ia 

We want to apply this result to the functions S x u x , d in Definition 12. II By (|2"3"|) we 
obtain 

Corollary 3.2. a) Let (o~,p,q) verifying 
n 1 n 

a = 1 h -, 2<p<oo, 2<<7<oo 

2 p q 

Then for (cr,p,q) as in (|37[) we have 

\\S x Vu\\ LPLq < \ 1 -'-°dl- 6 \\u\\ x ..9 
If additionally 9 > | f/ien 

||5aV«||lpl« < A 1 — 
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b) If instead 

2 n- 1 n- 1 

- H > 

P 9 " 2 



The interesting triplets of indices for (a, p, g) in 4 + 1 dimensions are 

(0,oo, 2) (energy) H ) H)(Strichaxtz) (-|, 2, 6)(Pecher) 

2 J J 6 

In addition, we can also use the index q = oo. Thus we obtain the triplets 

3 

(-2,oo,oo), (--,2,oo) 

For the case when < |, we rely on the additional triplets 

(-^2,3), (1,4,2) 

For convenience we summarize the bounds we need for AT^'^: 
Corollary 3.3. For < 9 < 1 we have 

x^WSxVuhoo^ + x s -i\\s x vu\\ L 2 L oo + x s - 3 \\s x vu\\ L oo <d*- fl |M|* 

A-^HSaVuII^ls + A-^ISaVuIU^ < di-°||«|| A .. 9 

A ,d 

The reason we include the gradient is to have also bounds for Ut ■ Because of the 
frequency localization, if we drop the gradient the same bounds hold with one less 
power of A. 

In our estimates later on we also need to work with X s ' b functions which are 
concentrated into a smaller modulation range. For this we introduce the additional 
norm 

{d d ~\ 

V ||u/i|||. s ,e ; M = Vm 
A , h 1 
h=l h=l ) 

If d = A we simply write X^' e . A simple argument leads to 

{00 00 ^ 

^2 WSxuxW^s.e; u = ^2s x u x > 
A=l A=l J 

We also have 

Corollary 3.4. a) Assume that 9 > \. Then 

A-^ISaVuIUoo^ +\ s -i\\S^Ju\\ L ^ +A s - 3 ||5aV u || L oo < 

A , < d 

b) Assume that 9 < |. Then 

A'-^ISaVuIUco^ + \ s -?\\S x \7u\\ L2L ^ + \ s - 3 \\S x \7u\\ L ~ < d?- e \\u\\ x*,e 
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In preparation for proving bilinear estimates for the X s,e spaces we first investi- 
gate which multiplications leave the Xt'j space unchanged. For this we define the 
algebras Md, M<_d with the norms 

ll/lk, = H/IU- + cT 1 ||/ f || L ~ + <H||/t|| £ »z- + d-i\\f tt \W L °° 

\\f\\M <d = \\f\\M d +^\\f\\ L ^ 

Then we have the multiplicative properties 

Lemma 3.5. Assume that f is localized at frequency d < A. Then we have 

(39) WfSxu\\x«,° < \\f\\M d \\u\\^,s 
respectively 

WfSxu\\ x ^ >< \\f\\M <d \\ u \\xV e > 9 <\ 

(40) x - d A<d 1 

\\SSxu\\ ± , A <d«-hf\\ M< AA\xv» ' e> \ 



The proof is straightforward, using Leibnitz's rule and the energy estimate (f26j) . 
To bound functions in the Md, respectively M < d norms we use Corollary 13.41 with 
d = A to obtain: 

Lemma 3.6. a) Assume that 9 > i. Then 

H^A^H a/ <a < A 2 - S ||w|| ^, fl H^auIIma < max{l, X 2 - s }\\u\\ Xs ,s 

\\S <x u\\m <x < max{\^ , \ 2 - s }\\u\\ X e,e 
b) Assume that 9 < i. Then 

\\Sxu\\ M<x < \^- e - s \\u\\ x s,e \\S <X u\\ M , < max{l, \i~ e ~ s }\\u\\ Xs ,e 
I|5'<A'"||m<a < max{A2, \^- e - s }\\u\\ X s,e 

Using the above property we prove the algebra property for the space X. 
Proposition 3.7. Assume that s > 2 and | < 6 < s — | . Then X s ' is an algebra. 

Proof. Let u, v € X s ' 6 . For both we consider the decomposition in Definition 12. 1[ 

oo A oo A 

u=^^s x ux td , v = y^ y^ S\v\,d, 

A=l d=l A=l d=l 

For the terms in the decomposition we use the XV d norms > as allowed by Corol- 
lary [231 We denote 

A d 
"A = ^WA,d, "A,<rf = y^^A./t 

Then we write 

oo oo oo oo 

UV = ^^(w) = S ^( S M u XiSx 2 vx 2 ) 

M=l /i=lAi = lA 2 = l 
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There are two cases when the above summand is nonzero, namely if Ai « A2 > /i 
and if max{Ai, A2} ~ /i. We consider them separately. 

Case 1, Ai, A2 ~ A > //. In this case the summability with respect to A is 
trivial, so it suffices to look at the product S\u\S\v\ for fixed A. This is localized 
at frequency < A. Combining the L°°L 2 and the L 2 L°° bounds in Corollary |3.4l we 
obtain 

(41) \\S x u x S x v x \\ L 2 + X-^dtiSxuxSxv^Wv < A- 2s+ t|| WA ||^a, e |H|^, e 
Using the equation we can also bound the second time derivative, 

(42) \- 2 \\d 2 (Sxu x S x v x )\\ L 2 < A- 2s +§|K||^, 9 |Hlx^ 
The three bounds above allow us to estimate for fj, < A 

\\S x u x Sxvx\\ x ..* < ^ s+e - 2 X- 2s+ i \\ux\\x-.» hxWx'S 

This suffices provided that 9 < s — §, which is insured by our hypothesis. 

Case 2. Here we consider products of the form S^v^SxUx where fi <C A. Then 
the product is localized at frequency A. The summation with respect to A is trivial, 
but not the one with respect to fi. We write 

A 

S^v^Sxux = S^v^Sxux^n + S^v^Sxux.d 

d—^i 

Using Lemma 13.51 and Lemma 13.61 we obtain 

A 

WSxUxS^V^y^e < \\SuVuSxUx,<u\\%-s,B + \\S l _ l V l _ l SxUx,d\\%s,e 

X X.Ji A X,d 

d—fi 

A 

£ ^^WS^V^Wl^Jux.^W^s.o + II^^Hl^ \\ u \d\\ 2 xs,e 



d—p. 



d=l X ' d 
< M 3 + 2e - 2s ||^|||, e ^|| UA , rf ||| s , e 



d=l 

The summation with respect to [i is trivial since 9 < s 



□ 



We next prove (|T3j) . 

Proposition 3.8. Assume that s > 2 and \ < 9 < s — | . TTiera we Ziave i/ie 
multiplicative estimate 

~£sfi _ j^-s-i,e-i j^s-1,0-1 

Proof. By duality this reduces to the multiplicative estimate 

xs ,0 . (x i- s ,i-e + L 2 H 2-s- e) c x i- s ,i-e + L 2 H 2-s-e 

Since s > 2 we have the fixed time multiplication 
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which implies the space-time bound 

Due to the energy estimate for X s ' it remains to show that 

](s,6 . 8,1-9 (- j£-l-s,l-0 i ^2 j^2-s-0 

We consider a product SyaxS^v^ which we decompose as in the previous proof. 
Because of the lack of symmetry we now need to consider three cases. 

Case 1. Here we estimate S^(S\u\S\v\) where /i < A. By Corollary 13.41 we 
obtain 

\\S\U\S\V X \\ L2L i < \\S\U\\\ L 4 L z\\S\Vx\\n L 3 < ^^iWuxW^eWvxW^s.i-e 

Using then Sobolev embeddings we obtain 

\\S^S x u x Sxvx)\\l^— o <^-'- e \ e -Hu x \\^B\\vx\\jtx-.,x-, 

Case 2. Here we bound S^u^Sxvx, The product is localized at frequency 

A, and the analysis is almost identical to Case 2 in Proposition 13. 71 

Case 3. Here we bound SxUxS^v^, (i <C A. The same argument applies, the 
only difference is that we gain some extra /i/A factors. 

□ 

We continue with the Moser estimates in (|T2]) . which follow from 

Proposition 3.9. Assume that s > 2 and i < 9 < s — | . Let T be a smooth 
function. Then 

l|r(«)|U..» <c(H L ~)(i + H&.. fl ) 

Proof. We write 

T(u) - r(v) = (u - v)f(u,v) 

and 

f(u,v) - f(x,y) =(u- x)g 1 (u,v,x,y) + (v - y)g 2 (u,v,x,y) 
where /, g\ and g 2 are smooth functions. Then we have 

oo 

r(u) =r(«i) + r(«<2A) - r( u < A ) 

A=l 

oo 

=r(ui) + ^ u 2\f{u<2\, U<x) 
A=l 

oo A 

=T(ui) + ^2u 2 x[f{u<2,u 1 ) + ^(/(u<2 M ,'u< M ) - /(n< M , u<^ 2 ))} 

A=l M=2 
oo A 

=r(ltl) + ^U2a[/( u <2,Ui) + ^( U2 M.9l( U <2M> U <M' M <A'/2) 
A=l M=2 

+ Ufi 52("<2 M , u< M , u<n/u))] 
Hence we need to bound expressions of the form 

SxuxSftVfthiS^w), i-i < A 
There are two different cases to consider: 



16 



DAN-ANDREI GEBA AND DANIEL TATARU 



Case 1. fi w A. Then the product has the form 

S\u\ S\v\ h(S<:\w) 

The first product is localized at frequency A and can be estimated as in (|41"|) . (|42l) . 
For the nonlinear expression we use Lemma 13.61 to obtain 

II<S<aHIm a < \\w\\ X s,e 

On one hand by the chain rule we obtain 

(43) \\h(S <x w)\\ M , < C(|HU-)(1 + \\wf xs , e ) 

On the other hand because of the frequency localization we also have the improved 
high frequency bound 

(44) \\S^h(S <x w)\\ M ^<C(\\w\\ L ~)(-) (1 + IM&...), M»A 

Taking this into account and repeatedly using Leibnitz's rule we get 
\\S\u\ S\v x h(S < \w)\\ xs , e 

oo 

~ £ II^C^AWA S\V\ h(S < \w))\\ 2 j i s,e 

< ^2 \\S\uxS\v x h(S < xw)\\ 2 i ^ e + ^2 \\S\u\S\v\S M h(S < \w))\\j l ., l g 

<c(||HU~) f E m 2s+28 - 4 a- 4s+7 + £ x w+3 - 2s (±Y) KII^..|M&..(i + | 

< C(\\w\\ L ~) A 2e+3 - 2s 11^11^11^11^(1 + ||«,|&.,.) 

This is trivially summable with respect to A. 

Case 2. /i « A. Then the product has the form 

S\u x S^v^ h(S <fl w) = 
Sxux, <fJ , S^iv S <fl h(S <fl w) + £ Sxu\^ < dS f _ l v t _ l S d h(S <ll w) 

/i<d<A 

+ £ Sxux^S f _ l v^S <d h(S <fl w) + SxuxS f _ l v f _ l S\h(S <ll w) 
+ £ 5 a ma S^h(S <tl w) = fi + h + f3 + fi + h 

^»A 

For fi we use Lemma 13751 Lemma T3. 61 and (|43|) to obtain 

ll/i|lx s - 8 ~ \\Sxux,<^S^v^\\ xS ,e \\h(S <fl w)\\ Mtl 

A,U A,I£ 

< / _ 5|| UAi<At || . e ||S>Jm<J w)||m„ 



< c(||HU«V + *- fl IK< M || xi>° \M\ x-A 1 + 

The summation with respect to fi is trivial and the square summability with respect 
to A is inherited from the first factor. 

For /2 we apply the same argument. There is a loss of a small power of {d/ fif 
from the first product, but this is compensated by the gain of arbitrary powers of 
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n/d due to |44|) . The same works for f% but there is no (d/n) loss. In the case of 
/4 we need to worry about the A summability, but the (fi/X) gain in (|44|) settles 
this. Finally, for /s there is a (/j,/v) n gain which cancels again all the losses. 
Summing up the pieces we obtain 

\\S u (S\u\ S^v^ h(S <f _ l w))\\ Xs ,e 

< C(||H|L~)^ +e - 2 A- 2s+ 5(^) w || WA ||^. e ||z;J|^,(l + \\ W \\ X .,,) 
for ^ <C A, 

||S , 1/ (5 A «A J S F M w /*M'S<M» y ))IU-* ^ C7(ll«'IU~)/* fl+f ~ S |l«A||jc ; .»ll«'Mlljc-.«( 1 + 

for ss A, respectively 

II^C^A^A /zC^^-^))!!^.^ < C7(||z i ;|| i _) M e +t- s (^) JV ||^ A || ||^ M || ^.^(l-t-II^H^.^) 

V AM 

for !/ > A. 

This concludes the proof of the proposition. 

□ 

Finally, we consider the bilinear estimate in (|14[) , which follows from the next 
Proposition. Its proof cannot be completed using the type of arguments we have 
employed so far. Instead, we contend ourselves with reducing it to the trilinear 
estimate in (|50|). to the proof of which we devote the rest of the paper. 

Proposition 3.10. Assume that s > I and § < 6 < s — § . Then we have the 
multiplicative estimate 

(45) ||VmVu||x«-m-i < IMIx-«IMIx 3 ' 8 

We begin our analysis with a simple observation, namely that 

Lemma 3.11. If u € X s ' 8 then Vu € X s ^ lfi where 

Xs-i,b = X s-i,e + ( L 2 H *+e-i n H i H s+s-2y 

Proof. We first consider spatial derivatives, for which we prove the better bound 

\\W x u\\ Xs -i.e < \\u\\ X :e 

By Definition 12.11 and Corollary [23] it suffices to show that for functions v localized 
at frequency A we have 

l|V x w|| xr i,e < \\v\\ xV » 

A , a A , a 

But this follows from the straightforward commutator bound 

(46) IIPff <vx ,VMU a <A||V«|| L . 

Here we recall that g 00 = 1, therefore every term in the commutator has at least 
one spatial derivative. 

Next we consider time derivatives, where it suffices to show that for functions 
v localized at frequency A we can write v — v\ + v 2 where v\ , 1)2 have the same 
frequency localization and 

/a y~ e 

(47) \\d t vi\\ x s-i,e + f - J ||^U2||( i 2 ff .+«-inH 1 H*+< , -») ^ IMI^m 
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Roughly speaking Vi accounts for the low modulation (< A) part of v while v 2 
accounts for the high modulation part. We define v 2 as 

This satisfies the bound 

\\w 2 v 2 \\ L2 < \\a g<v _v\\ L2 

which implies both the v 2 bound in (I47| and an H 2 bound for v\ which gives the 
correct L? bound for d t v\, 

A s_1 d 9 ||a t ui|| L 2 + ( - J \\dtV2\\( L 2 Hs +e-i nH i Hs +e-2 ) < \\v\\x s x - e d 
It remains to estimate □„ ,-dtV-i . We have 

Wn^dtv^w < Hp^.ftKll^ + \\d t n g<y/x vi\\ L * 

For the first term we use again (j4l)|) . For the second we compute 

D 9<v ^i = i- a 9< ^ + Ax, t )(Ax, t )- 1 n 9<vx w 

Since the difference \3 g< ^- — A x j contains no second order time derivatives this 
yields the bound 

\\d t n g< ^ Vl \\ L2 < \\\n g< ^v\\ L * 

This allows us to conclude the proof of (|47j) and therefore the proof of the lemma. 

□ 

We now return to the estimate Using the duality in (f3"2")) , (|l5|) reduces to 



(48) 



uvwdxdt 



< 



IMIxs-i.slMlxs-i.elM!* 1 - 3 ' 1 -^ 2 ^- 3 -" 



We do a trilinear Littlewood-Paley decomposition. Due to symmetry, we need to 
consider two cases. 

Case 1. Here we consider high-high-low interactions and bound 



I = J S\u S\v S^w dxdt 1 /j < A 



We have 

1^1 ^ ll'5 , A'ti||L~L 2 ll' S 'Al'||L 2 L s ll 5, A' u 'IU 2 L 3 

which by the embeddings in Corollary 13. 21 give 

\I\ < A*- a '+ a /i' +fl -S||«|| jf ._ 1 ,,||«||^_ 1 ,.||ti;|| Jcl -.. 1 - <+£afl »-.-. 

This suffices since both the exponent of A and the sum of the two exponents are 
negative. 

Case 2. Here we consider high- low-high interactions and seek to bound 



I = J S\u S^v S\w dxdt, fi <C A 



As a first simplification we dispense with the auxiliary L 2 norms. Begin with 
\I\<\\S x u\\ L 2\\S^v\\ L -\\S x w\\ L 2 

;$ X s ~ 1 /j,i\\S\u\\ L 2 M «->||£ S _ M A^HSaHIl 2 

This allows us to dispense not only with the L 2 H S+B ~ 1 part of u, but also with its 
X\~>^ component. 
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If v € £ 2 iJ s + 6 '- 1 n H X H B+6 - 2 then we bound 

1^1 < \\S\u\\ L oa L 2\\S u v\\ L 2 L ^\\S\w\\ L 2 

<m 8 — 9 IMI*.-m MM*.-! X^WSxwh* 

Finally, if w G L 2 H 2 ~ s ~ e then we can also estimate 

\i\ < I|SaHU~lH|£m 1 'IU 2 l~II>S'aHU 2 

which suffices for both the L 2 H 2 ~ s ~ e and the X)J>u~ components of w. Hence 
we have reduced P8|) to the bound 

(49) 



1^1 ~ ||S'Au|| xS -i,«||S' M w|| X3 -i,e||S , Aw|| J( -i- s ,i-e /i < A 



Unfortunately we cannot fully prove this using Strichartz type estimates. How- 
ever, we can use scaling to simplify this further and reduce it to 



(50) 



S\u S u v S\wdxdt 



<\nn |M| x o,i||u|| s ilMI o,l M< A 



For now we show that (|50p implies (|4"9"| . The remaining sections of the paper are 
devoted to the proof of (|50|) . 

After cancelling the powers of the high frequency the estimate (|4"5|) follows after 
summation with respect to 1 < di , d-i , cZ 3 < fi from the bounds 



(51) 

if c?2 < d. 
(52) 



S^w S^u S\wdxdt 
, respectively 



-Saw S^u S\wdxdt 



^ m /" Cin^mid^allwllx?'" NLf.o Hlx?'° 



X 



^ ln M dmi n dmaa:||«|| x o,o ||u|| 6 „ \\w\\ x o,o 



A \,dl X 4 ' A *,d 3 



To reduce all these cases to ([50]) we use scaling combined with a time decompo- 
sition argument. Precisely, for 1 < d < A we consider a smooth partition of unity 
in time with respect to time intervals of length 

Then a simple commutation argument shows that we can localize the Xt' d norm 
to the d~ l time intervals while retaining square summability, 



(53) 



3 



3 II 2 

Xd u 



x: 



We use such time decompositions in order to carry out the following three re- 
duction steps: 

(i) Reduction to d m i n = 1. By ((53)) all three norms are square summable with 
respect to time intervals of length Hence it suffices to prove the bounds on 
drain time intervals. Rescaling such time intervals back to time 1 we arrive at the 
case d m i n — 1- The regularity of the coefficients improves after the rescaling, here 
and below. Also we note that by Duhamel's formula we can replace the factor 
corresponding to d m i n by a solution to the homogeneous equation. 

(ii) Reduction to d m id = 1- By lp)3")) the norms corresponding to d max and d m id 
are square summable with respect to time intervals of length d^ id . Hence it suffices 
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to prove the bounds on d~^ id time intervals. Rescaling such time intervals back to 
time 1 we arrive at the case d m id = 1. Again by Duhamel's formula we also replace 
the factor corresponding to d m id by a solution to the homogeneous equation. 

(iii) Here we are in the case where two of the factors are solutions for the homo- 
geneous equation. In the case of (|5ip the remaining factor is at high frequency A; 
then we use directly (|50| . 

In the case of (f52j) the remaining factor is at low frequency /i, so we need to 
prove that 



S\u S^v S\wdxdt 



< ln/z d^\u\\ x o,o\\v\\ x 5 i<t \\w\\ x «.o 



Partitioning the unit time into about d time intervals of length d 1 this would 
follow from 



Xd,S\u S^v S\wdxdt 



< In d 4 \\u\ 



\W\\ v o,a 

" A A,1 



Rescaling the small time intervals to unit size this becomes exactly (|5 



4. Half-waves and angular localization operators 
We write the symbol for D g , 

/,/.,-. -.ii r 2 ! r^,C ! 

in the form 

p(t,x,r,0 = (T + a + (t,x,0)(r + a-{t,x,0) 
This leads to a decomposition of solutions to the wave equation into two half-waves: 

Proposition 4.1. (Geba-Tataru |16j ) Let u be a solution to the inhomogeneous 
equation (|30[) for D g . Then there is a representation 

X7u = u + + u~~ 

where 

\\u+\\ L 2 + \\(D t + A+(t,x,D))u+\\ L 2 + \\u-\\ L 2 + \\(D t + A-(t,x,D))u-\\ L 2 

^ \\ u \\m + \\O g u\\ L 2 



As a consequence, in ([50)1 we are allowed to replace solutions to the O g equation 
by solutions to the D t + A + , respectively D t + A~ equation. We also denote 

\\u\\x ± = \\u\\ L 2 + \\(D t + A ± (t,x,D))u\\ L 2 



M\x ±i 



d*\\u\\ L 2 + d-^\\(D t + A ± {t,x, D))u\\ L 2 



In order to facilitate the use of microlocal analysis tools it is convenient to replace 
the symbols with mollified versions a<„ defined by 



a <n(t, x >0 = S <t ,(D x )a(t,x,£,) 



Given an angular scale a we consider the ± Hamilton flows for D t + Az a _ t . 



(54) 
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These are bilipschitz flows, homogeneous with respect to the £ variable. The an- 
gular scale is relevant in that the Hamilton flow for Dt + -<4< a -i serves as a good 
approximation to the Hamilton flow for D t + A up to an 0(a) angular difference. 

To characterize the higher regularity properties of these flows is convenient to 
introduce (see [4]) a metric g a in the phase space, defined by 

ds 2 = |£r 4 (£^) 2 + ir 4 a- 2 (£ A d£) 2 + a-^\- 2 (^dx) 2 + |r 2 a- 2 (£ A dx) 2 
Then as in [1] we obtain 

Lemma 4.2. The Hamilton flow maps (xf ,£j ) are g a -smooth canonical transfor- 
mations. 

Given a direction 9 e at time t = we introduce the size a sectors 

S a (8)={Z; Z(Z,8)<a} 

S a (6) = {£_; Ca < Z{£_,9) < 2Ca} 

where C is a fixed large constant. The images of K" x S a (9), respectively R™ x S a (9) 
along the Hamilton flow for D t + A^ a _ 1 are denoted by H^S a (9), respectively 
H±S a (9). 

Let £g = £g(x,i) be the Fourier variable which is defined by the D t + ^4< Q -i 
Hamilton flow with initial data £#(x, 0) = 6 (i.e. (x, t) = is the solution of 

the flow (|54|) with initial data ^ — £, for which xf (t) — x). This is well defined 
at least for a short time, precisely for as long as caustics do not occur. From 
Lemma |4~21 one also sees that £g is a (7 a -smooth function of x. 

We consider a maximal set O a of a-separated directions and a partition of unity 
at time 

9eo a 

consisting of 0- homogeneous symbols supported in S a (9) which are smooth on the 
corresponding scale. Transporting these symbols along the ± Hamilton flows by 

produces a time dependent partition of unity 

(55) f = E xf' a (t,x,0 

9eO a 

so that the support of Xe'° "(^i^iO is contained in H^S a (9). 

The regularity of these symbols is easily obtained from the transport equations 
(see again [3]): 

Lemma 4.3. The symbols x^' a (t, x, £) belong to the class S(l,g a j^. 

We use the above partition of unity in the phase space to produce a corresponding 
pseudodifferential partition of unity. Given a frequency A > cT 2 we define the 
symbols 

Xe,\(ti x ,0 = S < x/s(D x )xp a (.t,x,^)S x (0 



^Throughout this paper we will use the standard notation S(m,g), while in [4] we used for 
5(1, g) the shorter one: S(g). 
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These are used in order to split general frequency localized waves into square sum- 
mable superpositions of directionally localized waves, 

S\u = ^2 Xe,' X {t,x,D)S x u 
9<=o a 

This decomposition is closely related to a wave packet decomposition, see [TT], [T2] . 
[TB] . and [3]. The difference is that here we skip the spatial localization part since 
it brings no additional benefit. The above localization at spatial frequencies less 
than A/8 insures that the output of the operators Xg '\{t, x, D)S\ is still localized 
at frequency A. This localization is otherwise harmless: 

Lemma 4.4. The symbols Xq' x {Pi x iC) belong to the class S(l,g a ). In addition, 
we have similar bounds for the Poisson bracket 

(56) {r + a* Q _ 1 (t,a;,^),x^ Q (i,a;,^)} G S{l,g a ) 



Proof. The fact that Xe'x^^iO € S(l,g a ) is straightforward since the multiplier 
S<\/8 18 a mollifier on the A" 1 spatial scale, which is less that the spatial scale of 
the g a balls. 

Since Xe' a is transported along the 0< a -i (t, x, £) flow, the Poisson bracket is 
expressed in the form 

{a ± <a - 1 (t, x, OMO}xt£(t, x, o + s A (0[ff a ± o _ i , S <x/8 (D x )]xp a (t, x, o 

Here H a ± is the Hamiltonian operator associated to the Q< a _i (t, x, £) flow, ft is 

easy to see that the first term belongs to S(l,g a ), therefore it remains to consider 
the commutator term. We have 

The commutator of a scalar function g with S <x / S can be expressed as a rapidly 
convergent series of the form 



[g,S <x/8 ] = X- 1 ^S L < :{ /s VgS, 



<A/8 



3 



where the multipliers S^, s and S^.j have the same properties as S<\/8 an d decay 
rapidly with respect to j. Then the above commutator term is expressed as 

[H^ , S <X/S (D X )} = A" 1 £ S]j{ /8 {d.d^.A - dla% a _^) Sf x/8 

j 

At this stage the effect of the mollifiers is negligible and we can use the regularity 
properties of a ± and Xe' a to directly compute 

~sx(0[H a ± ,S <x/s (D x )}xp a (t,x,Z) G S(^-,g a ) 

<a 1 QL A 

□ 

To better understand the phase space localization provided by Xo' x consider 
some point (x ,t ) and the corresponding center direction £g(x ,t ). A spatial 
unit g a ball Bf(xo, to) centered at (xq, to) has dimensions^ a 2 x a" -1 with the long 



2 Here n stands for the space dimension 
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sides normal to C$( x 0i to)- Within the ball Bg(xo, to), Xe'x&x localizes frequencies 
to a sector of angle a centered at £g(xo,to). Thus the frequencies are localized to 
a radial rectangle centered at A£# (xo,to) of size A x (aA)" -1 . In this picture, angle 
a wave packets correspond to a spatial localization on the scale of the above ball 
Bg(xo,to), constructed along a fixed ray of the Hamilton flow. 

The <7„ metric restricted to frequency A is slowly varying and temperate at fre- 
quencies □ A > oT 2 , and in our analysis we will always be above this threshold. 
Hence there is a good pseudodifferential calculus for operators with S(l,g a ) sym- 
bols. The semiclassical parameter h = h{a, A) in the S(l,g a ) calculus at frequency 
A is given by 

h(a,X) = (a 2 A)" 1 
The S(l,g a ) symbols at frequency A satisfy the bounds 

(57) \(^d x r(^ Ad x fd^(^rq(t,x,0\ <a- 2 ^(a\)- v 

Due to the L 2 in time regularity of the second order derivatives of the coefficients 
we also introduce the space of symbols L 2 S(l,g a ) which at frequency A satisfy 

(58) \(^d x r(^Ad x f(f^d e rq(t,x,0\ <or 2 °-W{a\T v f{t) 

for some /eL 2 . In all the operators we consider here, the function / is the same: 

(59) f(t) = M(\\V 2 g(t)\\ L ~) 

In some of our estimates we need to deal with two distinct scales at a given fre- 
quency A, namely the angular scale a and the A^ scale at which the coefficients 
are truncated. Correspondingly we introduce additional symbol classes C%S(l,g a ) 
of symbols q localized at frequency A which satisfy the S(l,g a ) bounds (fS"7|) for 
c + |/3 1 < k, respectively the weaker estimate 

(60) \(^d x n^Ad x fd^(^r q (t,x,0\ < (a- 2 °-W+a- k \^)(a\T» 

for <t+ |/3 1 > k. There is still a calculus for such symbols, since the above bounds are 
stronger than the S(l,g i ) bounds. The related classes of symbols L 2 C^S(l,g a ) 
are defined in a manner which is similar to (|58[) . 

Using the calculus for the above symbol classes one can prove that the partition 
of unity in (|55p yields an almost orthogonal decomposition of functions, namely 

Proposition 4.5. Fix a frequency A and let a > A~5 . Then for each function u 
which is localized at frequency A we have 

(61) E \\xix(^D)uf x± ^\\uf x± 
oeo a 

Proof. We only outline the proof, since this result is essentially contained in [15] . 
There are two bounds to prove. The first 

(62) £ \\ X ^(t,x,D)u\\ 2 L2 ^\\u\\ 2 L2 

0£O a 



J This corresponds to the classical wave packets which are localized on the scale of the uncer- 
tainty principle. Above this threshold we are dealing with generalized wave packets, which may 
have a more complex structure, see I16| and [2] 
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follows from the almost orthogonality of the operators Xg '\(t,x,D). This in turn 
is due to the almost disjoint support^ of Xe'\ an d to the 5(1, g a ) calculus. 
Consider now the second bound 

(63) II ^ + A± )xtx(t, x, D)uf L , « || (A + A^uWl* + 0(\\uf L2 ) 
eeo a 

We first establish it with replaced by A 1 , 

<A 2 

(64) Y, \\( D t + A ± ^f(t,x,D)u\\l 2 ^\\(D t + A ± )u\\l 2 +0(\\u\\l 2 ) 
eeo a 

Due to l|62p and the energy bound 

\\n\\l^<\\u\\h + \\u\\ L A\{Dt + A ± L )u\\ L . 

<A 2 

it suffices to prove the commutator estimate 

(65) \\l D t+ A± xi ,xtx(t^,D)Ml^ IMI£«l* 
eeo a < 2 

which we split into two components. 

For the low frequency part of the coefficients we use a second order commutator 

(66) Y \\l D * + K„-i>x*£b*> D )Ml> Z Nl£~i> 

eeo a 

For this it suffices to prove that 

(67) [D t + A± a _ 1 , xt'x(t, x, D)] e OPL 2 S(l, g a ) 

The summation with respect to 9 G O a follows by orthogonality since the symbols 
for the above commutators will retain the rapid decay away from the support of 
Xe\- Here it is important that (|59[) applies uniformly. 

Due to the Poisson bracket bound in (|56p it suffices to show that 

[^<a-> . xt'x(t, x, D)] + i{a± a _ x , xtf}^ D ) e OPL 2 S{\, g a ) 

Due to the frequency localization of Xq'x ■> only the values of a ± (t, x, £) in the region 
|£| i=a A can affect the above operator. At this point it is no longer important that 
a^ Q _! and Xe'x are related. We consider a rapidly convergent spherical harmonics 
expansion of a , 

a ± (t,x,0=Y b ^ x ^(0 

3 

where bj have the same regularity as the coefficients g^ while <f>j (£ ) are homogeneous 
of order 1. ft suffices to consider a single term b(t,x)cj>(£) in this expansion and 
show that 

(68) [b^-.^^tx^^^+^a-^xtxH^^D) e OPL 2 S(l,g a ) 

To see this we consider the commutators with b and with </>. The commutator term 
with b has the form 

C b = ([&< : (t,x, D)]+i{b < 

a" 1 ) Xg\ 

l }(t,x,D))^(D) 



modulo tails which are rapidly decreasing on the g a scale 
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Since Sg& <ot -i € L 2 5(l,5 a ), dfx^f € S^A- 2 , ffa ) and e S(X,g a ), the S(g ) 
calculus at frequency A yields the better result Cf, S OPL 2 S(a~ 2 \~ 1 , g a ), which is 
tight only when a = A~ 2 . 

The commutator term with <f> has the form 

= b <a -i MD), X ^(t, x, £>)] + i{cf>, xt;x}(t, x, £>)) 

The b <a -i factor belongs to S(l,g a ) and can be neglected. The argument for the 
remaining part is somewhat more delicate since it hinges on the homogeneity of <j>. 
With b = 1 denote by £ the input frequency for and by r\ the output frequency. 
Due to the homogeneity of (f> we have the representation 

(69) 4>( V ) - HO = (r? - OV0(O + A (£ — t?)) 2 

where -0 is a smooth and homogeneous of order —3 matrix valued function. For 
|£|, \r)\ « A we can separate variables in ^ and express it as a rapidly convergent 
series 

This gives a representation for of the form 

C^A^^^WC^A^x^Ct.x.D^P) 
i 

Since Xg '\( X >D) G S(l,g a ) we obtain (£ A d x ) 2 Xg'\ S 5(A 2 a -2 , g Q ) which shows 
that € OPS(a~ 2 X^ 1 ,g a )- This concludes the proof of ([68]) and thus the proof 
of (EH). 

For the intermediate frequency part of the coefficients we have a first order 
commutator estimate 

(to) Yl Ht A± , ^,xtt(t,x,D)Hl* < 

eeo a 

Together with |66|) this implies (|65]). 

This follows from first order commutator estimate 

(71) [A± ljX ±n^^)]eOPL 2 Ci5(l l!?a ) 

a 1 <-<A 2 

Indeed, for a scalar function b we can estimate 

°" 9 |l 6 a-i<.<A*U^- +a! " 1 |l^a-i<.<A*ll^- ~ 

Applied to the the symbol a ± as a function of x this shows that 

a* , EL 2 ClS(a 2 \,g a ) 

Since x^'x e S(l 7 g a ), the estimate (|7Tj) follows by pdo calculus. The square 
summability with respect to 6 is again due to the almost disjoint supports of the 
symbols Xg' a - 

It remains to pass from (|64|) to (|63|) . Due to the energy bound 
IM| 2 l~ l2 < ||u||i a + \\u\\ L 4(D t + A±^)u\\ L * 
this is a consequence of the estimate 

11^ i«IU»<HU-L» 

>A2 
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applied to both u and Xg'\(t> x > D)u. Using the spherical harmonics decomposition 
of the symbols a ± as above this reduces to the straightforward bound 

\\b >xi u\\ L 2 < X-^bh^WuW^ 

□ 

The frequency localization in Xg \ contributes to improved Strichartz type esti- 
mates above the critical range of exponents. Begin for instance with the endpoint 
L 2 L 6 Strichartz estimate 

(72) \\xt£(t,x,D)u\\ L z L e<\i\\u\\ x ± 

Here the angular frequency localization plays no role. However, suppose we want 
to use Bernstein's inequality to replace this by an L 2 L°° estimate. Modulo rapidly 
decaying tails, within each spatial g a ball Bg(xo,to) the function Xe'\(t> x > D)u is 
frequency localized in a dyadic sector section of size A x (aX) 3 . Then the constant 
in Bernstein's inequality is 

[A x (aA) 3 ]s = \ia^ 
Hence we obtain the better L 2 L°° bound 

(73) \\ x $;Z(t,x,D)u\\ L * L oo <a'A*||u||x ± , a > X~i 

A simpler related uniform bound is derived directly from the energy estimates, 

(74) \\ X ^(t,x,D)u\\ L oo< a ix 2 \\u\\x ± , a > A"* 

A similar bound holds for the right hand side of the Xg'\(t> x i D)u equation. Indeed, 
for u £ X± we can write 

(A + A ± )x^f(t, x, D)u = (D t + A±^)xt£(t,x,D)u + A± xh xtx(t, x > D > 

The first term belongs to L 2 and has a similar frequency localization as Xg'xi^i x i D)u. 
The second is estimated directly using ([73")) . This yields 

(75) \\(D t +A ± ) X ^(t,x,D)u\\ L 2 L ~<aix 2 \\u\\ x± , a > A"' 

Another way of taking advantage of the angular localization is in corresponding 
bounds for derivatives. Consider the differentiation operators A D whose symbol 
vanishes in the £g direction. Then in the support of Xg '\ these symbols have size 
aX. Hence from (|72p we also obtain 

(76) \m AD) X g h X(t,x,D)u\\ L , L e < (aX)xi\\u\\ x± 

We can argue in the same way for the energy estimates or for the L 2 L°° bound in 
(|73| . For convenience we collect several such bounds in a single norm, 

\\v\\ x x. a ,e = \\v\\ x± + \\v\\ l °°l 2 + *~^NU 2 £ 6 + a~^X~^\\v\\ L 2 La o + a^A~ 2 ||v|| L oo 
+ a-iX- 2 \\(D t + A^vWl*^ + (aA)- 1 ^ A D)v\\ L *o L z 
+ (aX)~ 1 (X~^\\{Cg A D)v\\v L * +a-U-f ||fe Q A D)v\\w) 

and use it to state a corresponding version of (|6ip , 

(77) J2 \\xtx(t, x ,m\W,*™\\Sxuf x± 
eeo a 
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We want to replace the partition of unity in (|55j) first with a bilinear one and 
next with a trilinear one. Given two frequencies fi < A, we denote a M = 
and introduce a corresponding bilinear partition of unity which is useful when 
estimating the frequency /i output of the product of two frequency A waves. The 
main contribution corresponds to opposite frequencies £ and 77, therefore we organize 
the following decomposition based on the dyadic angle a M < a < 1 between £ and 
—77. Precisely, by superimposing the a angular decompositions for a in the above 
range we obtain 

S\(£,)sx{v) = 

\6 1 +e 2 \<2Ca l , |0 3 +04|<4Ca: M 

E E xt$r (*> 0xl;7 (*, *7)4'?" (*- Oxg?* (*> »?) 

1 Ca<|e 1 +e 2 |<2Ca |e 3 +e 4 |<4Ca 

+ E E E xf; a x(t,x,OxJS(^ x ^)xe' 2 x(^ x 'OxJ; 2 x(t,x,v) 

To shorten this expression we redenote factors and harmlessly simplify the summa- 
tion notations to 

1 

(78) 1= *M*(^eW-#(*.*>'7)+ E E <nM,6<^A(M,?7) 

eeo QfI eeo a 

where the tilde in 4>^'" indicates an 0(Ca) angular separation from 9. The symbols 
(j>g'\ , respectively 4>^ retain the same properties as Xg '\i namely 

(79) 4;«eS(i, 9a ), {T + a± a _ 1 (t,x,Z),4;™(t,x,0}eS(i,g a ) 

and the same for 4>g'x ■ I n particular the counterpart of (|77|) is still valid, 

(80) E (*> x ' D >\\xf + Ui'x^ *> D H 2 x^ « \\Sxu\\ x± 
eeo Q ± ± 

Finally, we arrive at the main trilinear symbol decomposition. Its aim is to 
achieve a simultaneous angular decomposition in trilinear expressions of the form 

J uvwdxdt 

We denote the three corresponding frequencies by £, 77 and £. We assume that each 
of the factors has a dyadic frequency localization, 

|£|«M«A, Klra/i, l«/i<A 

If the trilinear decomposition were translation invariant then only its structure on 
the diagonal £ + r/ + £ = 0is relevant. However, in our case we are working with 
variable coefficient operators therefore a neighborhood of the diagonal is relevant. 
The size of this neighborhood is determined by the spatial regularity of the symbols 
via the uncertainty principle. 

Corresponding to the first term in (|78[) we consider a decomposition in £ with 
respect to the dyadic angle between £ and 9, 

mo = (*. x, + e x > 
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To understand the £ decomposition corresponding to the second term in (|78|) we 
first identify the location of the diagonal £ + i] + £ = 0. Given the above dyadic 
localization of £, r\ and C, if the angle between £ and — rj is of order a, then the 
angle between £ and ±£ must be of order aA/i -1 which is larger than a. Thus 
the interesting angular separation threshold for £ is aA^ _1 . It would appear that 
there are two cases to consider, namely when the angle between £ and C is small, 
and when the angle between — £ and £ is small. However, due to our choice of 
the ± signs corresponding to £, 77 and C, the latter case leads to nonresonant wave 
interactions and loses its relevance. Hence, the significant dyadic parameter here is 
the angle between £ and £, and the £ decomposition has the form 

/3>a/i- 1 A 

Then the full trilinear decomposition has the form 

eeo Qt , 

+ £ E fc a M>o 

(8i) +E E <t>t£( t > x >04-&( t > x >y)$t£ 1 ' x ( t > x >0 



E E 



+ E E ^a^'^'^^-m E ^V/f^'O 

In the above sum the first three terms are the main ones, as they account for the 
behavior near the diagonal. The remaining terms have off diagonal support, and 
their contribution to trilinear forms as above is negligible. 

5. Proof of the trilinear estimate (1501) 



As noted in the previous section, we can replace the spaces X^f d in (|50p with the 
X± spaces. Hence we restate (|50|) in the form 

Proposition 5.1. For any choice of the ± signs and 1 < d < < A «ie have 



(82) 



S\uS\v Snwdxdt 



< ln/x • \\Sxu\\x ± \\S\v\\x ±td \\S t j.w\\x ± 



Proof. We begin with several simple observations. First, by localizing to a fixed 
smaller space-time scale and rescaling back to unit scale we can insure that the 
coefficients g n vary slowly inside a unit cube, 

\Vx,tg ij \ « 1 

This in turn insures that the Fourier variable does not vary much along the Hamilton 
flow, 

\$ - 0\ « 1 

We can also localize all factors in frequency to angular regions of small size, say 
< ^j. The corresponding localization multipliers are easily seen to be bounded in 
X± and X± d- 
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If the first two ± signs are identical then the product S\u S\v is concentrated at 
a time frequency of the order of A which makes it almost orthogonal to S^w, hence 
the estimate above is much easier. Therefore without any restriction in generality 
we fix the first sign to + and the second one to — . Even though the problem is not 
symmetric with respect to the first two factors, the sign in the third factor plays 
no role whatsoever, so we fix it to +. We denote 

a(t,x,X) = a + (t,x,£) 

Then 

We note that, for the purpose of the above estimates, in the definition of X± at 
frequency A we can replace the symbols a(x,^) with their regularized versions, 
namely a <x i (x,£,). 

To keep the number of parameters small we first present the argument in the 
case when d = 1. Once this is done, we show what changes are necessary for d > 1. 

Case 1: d = 1. Corresponding to the trilinear symbol decomposition ([81]) of 
the identity we consider the corresponding pseudodifferential decomposition of the 
trilinear expression in (|82p . The we estimate each of the five terms. We remark 
that, since S\u, S\v and S^w are frequency localized in a small angle, so are all 
the factors in (f8"2|) . 

Case 1, term I: 

1 = I <f>t,f*{t,x,D)Sxucf>Z8 a £(t,x,D)Sxv 4>+'^(t,x,D)S fl wdxdt 

We use the energy estimate for the first two factors and the L 2 L°° bound for the 
third to obtain 

1^1 < II^a 1 " (^^ i^)^^!!^^^^,*) 1101^ C^' i^)^^!! ^A.^^.e H^^ 1 " ^, z?)^^!!^.^^ 

The summation with respect to 9 is straightforward due to (|80f . 
Case 1, term II: This is the most difficult term, 

H= I <t>tx' l (t,x,D)Sxu<j)Ze a x{t,x,D)S x v £ ^(t,x, D)S^w dxdt 

eeo a „ a>a„ 



The summation with respect to 9 is easily done using (|80p . Hence, in what 
follows, we fix 9 and redenote 

u e = ( l ) t.x"( t > x ' D ) s * u > v e = <t>Z'£x(t,x,D)Sxv, w% = 4>g , ^(t,x,D)S fl w 

The factors ue and vq are frequency localized in small angles around 9, respectively 
— 9; Wg has a similar localization around ±9 provided that a<l. 

We denote by x, £) the linearization of x(t,x, £) with respect to £ 

around £ = a:). Since a i(t, £) is a homogeneous symbol of order 1, we 

have 

a <j(j i (t, x, £) = £<%a <M i (*, x, ^ ) 



Consider now the difference 



a i — a i 

<M 5 <M' 



It vanishes of second order on the half line Due to the uniform (nonradial) 

convexity of the characteristic cone {r + a i (t, x, £) = 0}, it follows that e is 
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nonzero when £ is not collinear with £ e M . Precisely, we can estimate it in terms of 
the angle as 

e(t,x,£)«|£||Z(£,^)| 2 
In particular in the support of the symbol (j)g the above angle has size a and the 
frequency has size fi. Henc^l 

e(t,x,() ~ a 2 fi, (t, x, C) € supp 

Here it may help to think of the constant coefficient case where £g " = 8, while 
a — a = |£| — $6*. We introduce a local inverse for e(i, x, £) in the support of <^J~ 
namely 

The cutoff symbol <^J~^* is similar to 4%'^ but has a slightly larger support and 
equals 1 in a neighbourhood of the support of 

As defined, the operator L(t,x,D) is not localized at frequency fi. To remedy 
this we truncate its output in frequency and set 

L = S^{D)L{t,x,D) 

The properties of the operator L are summarized in the following 

Lemma 5.2. The operator L satisfies the following estimates: 
a) fixed time L p mapping properties: 

\\L\\lp^lp < a -2 /* -1 , l<P<oo 

b) fixed time approximate inverse of A(t, x, D) — A(t, x, D): 

\\(A(t,x,D) - A(t,x,D))L - <j>(t,x,D)\\ L p-> L p < fi~i + a" 2 //" 1 , 1 < p < oo 

c) space-time X + mapping properties: 

\\L\\x + ^x+ < or" 1 ^ 1 

Proof. We first compute the regularity of the symbol e(t, x, Q within the support 
of I. With respect to £ this is smooth and homogeneous, therefore we only have to 
keep track of the order of vanishing when £ is in the £ e M direction. With respect 
to x there is the dependence coming from the symbol a, as well as the dependence 
due to the £^ M direction occuring in the linearization. Since a is Lipschitz in x and 
is Lipschitz in x and smooth on the scale, within the support of I we obtain 

(83) eeC>V,S«) 

Combining this with the regularity of the symbol 4>$ € S(l, g a ) we obtain the 
symbol regularity for I, 

(84) l£ClS{{c?n)-\g a ) 

To prove part (a) of the Lemma we observe that for fixed (t, x) the symbol 
l(t, x, £) is a smooth bump function of size (a 2 ^) -1 in a rectangle of size fix [afi) 11 ^ 1 



here we switch to the letter f for the frequency, as the following analysis refers to the region 
at low frequency /i corresponding to the last factor w in the trilinear form. 
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oriented in the £g M direction. This implies that its kernel K (t, x, y) is bounded by 
(a 2 /i) _1 times an integrable bump function on the dual scale, 

\K(t, x,y)\< (a 2 M)-V(aM) n - X (l + (*> " + H^* (*. *) A (a - y)])"" 

This bound is symmetric; indeed, since ^^(t^x) is Lipschitz in x we can replace it 
by £g M (t, y) in the above bound. Thus integrating we have 

sup / \K(t,x,y)\dy < (a 8 /*) -1 , sup / |A"(t,a,y)|da; < (a 2 //)- 1 

x J y J 

The L p bounds for L(t, x, D) and also for L immediately follow. 

For later use in the proof we observe that within the support of I we have 

l£ M (t,aOA£| <afi 

Then the same argument as above yields the additional bounds 

(85) \\(C'(t,x)ADfLu\\ LP < (a/i^aV)- 1 

For part (b) we write 

(A(t, x, D) - A(t, x, D))L - 4>(t, x, D) = R x (t, x, D) + R 2 (t, x, D) 

where 

i?i(t, x, D) = E(t, x, D)§ li {D)L{t, x, D) - 0(t, x, D), 

respectively 

R 2 (t,x,D) = (A i(t,x,D)-A i{t,x,D))S^{D)L{t,x,D), 

The operator R\ is localized at frequency fi. The principal part cancels, and since 
e e Cj l S{a 2 ^, g a ) and I £ C^S((a 2 , g a ) by the pseudodifferential calculus it 
follows that 

R 1 (t,x,D)eC°S((a 2 ^- 1 ,g a ) 

In addition, the symbol of R\ decays rapidly away from the support of 4>q'^ ■ Hence 
we obtain the same kernel and LP bounds as in the case of L(t, x, D). 

Consider now the operator R 2 . Since a(t,x,C) is Lipschitz in x it follows that 
\a i (t, x, C)| < /x _T |C|- Expanding a i(t,x,() in a rapidly decreasing series 
of spherical harmonics with respect to £, we can separate variables and reduce the 
problem to the simpler case when a 1 (t, x, Q — b(t, x)c(() with \b\ < 2 and c is 
smooth and homogeneous of order 1. For the symbol c — c we use the representation 

c(C) - c(t, x, = V(C . K?" (*, x ) a C) 2 

where ^ is smooth in both arguments and homogeneous of order — 1 in £. Separating 
variables in ip we can assume without any restriction in generality that tf> depends 
only on £. Then after some simple commutations we obtain 

c(D) - g(t, £>) = (£" (t, as) A D) 2 ^(Z>) + 0(1) L p^ Lp 

To estimate this we use (|85|) . The factor ip(D)S^(D) yields an extra factor in 
the LP bounds, therefore we obtain 

\\R 2 (t,x,D)\\ LP ^ LP < fi-i 
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Finally we prove part (c). By (a), L is L 2 bounded with norm 0{a 2 fi 1 ), 
therefore it remains to prove the commutator estimate 

(86) \\[D t + A i (t,x,D),§rL(t,x,D)]\\ L ~ L ^ L *<ar 2 ir' 1 
This is a consequence of the operator bound 

[A + A i (t, x, D), S^L(t, x, D)] G L 2 C°5(a~ V\ 5a) 

2 " 

To prove it we use the pdo calculus to represent the commutator as a principal 
term plus a second order error, 

[A + A i (i, x, D),S^L{t, x, D)] = S^Qit, x, D) + R(t, x, D) 

where the principal part g has symbol 

q(t, x, £) = -i{r + a i (t, x, £),l(t, x, £)} 

The remainder R is localized at frequency /x. A direct computation, using (|84[) . 
shows that its symbol satisfies 

rei 2 C7°S(a-V\sa) 
It remains to consider the above Poisson bracket and prove that 

(87) qeL 2 C°S(a- 2 n-\g a ) 
For this we write q in the form 

T+,c* — 2 i -1 i -1 

l< l = -9e^^ e +92^ +q 3 e 

where 

?i = { T + a < M i' e }' 32(4, a:, f) = {t + a <Q ,-i > </>e,)T } 

respectively 

93(*>*>0 = {vi<.< M *'fe} 

Within the support of 4%'^ we know that e G C^S{a 2 n, g a ) is an elliptic symbol. 
Hence for the first term it suffices to show that qi <E C®S(a 2 /j,, g a ). Indeed, by 
definition q± is a homogeneous symbol of order 1 which is continuous in x and 
homogeneous in £. In addition, we know that e(t,x, C) vanishes of second order in 
£ at (t, x, £g M (a;,t)) which is also invariant with respect to the r + a i Hamilton 

flow. Then g must vanish of second order in £ at (£, x, £ e M (a;, i)). Arguing as in the 

case of e, this implies that within the support of ipg " we have qi € C®S(a 2 /i,g a ). 
As in ((67) ) we know that ao £ 5(1, g a ). Also we have a _ x< < i G L 2 C 2 S(a 2 fi, g a ) 

and </>+ p " G 5(1, 5 Q ) therefore q 3 G L 2 Ci5(l, g a ). 

This concludes the proof of ([57]) and therefore the proof of the lemma. 

□ 

To continue the estimate of term II in Case 1 we define the auxiliary trilinear 
form 



E(u,v,w) = J (A + A(t,x, D))uvwdxdt + / u(D t — A(t,x,—D))vwdxdt 
uv (A + x, D))wdxdt 
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With w = Lwg we write 

UgVgWgdxdt = 



ugvg ((A(t, x, D) - A(t, x, D))L - l)w^dxdt 
+ / (D t + A(t, X, D))tlg vgwdxdt 
+ / ug(Dt — A(t, x, —D))vg wdxdt 



ugvg(D t + A(t, x, D))wdxdt 
—E(ue,V0,w) 

We bound each term separately. For the first one we write 

(A(t,x,D) - A(t,x,D))L - 1 =(A(t,x,D) - A(t,x,D))L - 4>j'*(t,x,D) 

+ (}t£(t,x,D)-l) 

The contribution of the first line is estimated using Lemma (5TU(b) and (|75|) for wfi, 
u g vg[(A(t,x,D) - A(t,x,D))L - 4$£(t,x,D)]w%dxdt 

< (aT 2 /* -1 + ii~?)\\ue\\j J o OL ii\\ve\\L°°L*\\we IU 2 l°° 

< (a- 2 ^ 1 + Ll-?)\\ug\\ X+ \\vg\\x-\\w e *\\ L 2 L c X , 

< (cT 2 ^ 1 + H~^)a^ ^\\ug\\ X+ \\vg\\x_ \\Wg HjfM.c.8 

For the contribution of the second line we observe that 

&t£(t,x,D) - 1)< = (it£(t,x,D) - l^S^w 

where the symbols <^j~^* — 1 and ^e'^s^ have disjoint supports. Since they both 
belong to 5(1, g a ), this yields a gain of a factor (a 2 fi)^ N in (|75|) . with N arbitrarily 
large: 



E 



:(t,x,D)-l)w%\\l 2L „<»Ha^)- N \\ We <\\ 2 xfi ,, 



This is more than we need. 

For the second term in (|88p we use the L 2 bound for (D t + A)ug, the energy 
bound for vg and (|73p for w. This yields 



(D t + A(t, x, D))ug vgwdxdt 



< 



a 2 ^\\Ug\\ X+ \\Vg\\x.\\W B 



A" 



The third term is similar. 

For the fourth term in (|88|) we use the energy for the first two factors combined 
with Bernstein derived L 2 i°° bound for the third, 



ugvg (D t + A(t, x, D))wdxdt 



< \\u\\ x+ \\v\\ x _\\(D t +A(t,x,D))w\\ L 2 L ~ 

< (a 2 Liy 1 ( f i(an) 3 )^\\ug\\ x+ \\vg\\x_\\w e '\\ xlt , a ,e 
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It remains to prove the estimate for E. Observe that the time derivatives in E 
can be integrated out, producing contributions of the form 



(«» / 



uq vg wdx 



at the initial and the final time. These are estimated using energy bounds for the 
first two factors and the pointwise bound arising from Bernstein's inequality for the 
last factor, 

< (aVr'lKlk- < {a 2 »T\ l i{atf)Hwn xr ,<> = (^"M IKII^ 
This leaves us with a purely spatial trilinear form, 

J E (ug,Vg,w)dt 

where 

E a {u,v,w) = j A(t,x,D)uvw — uA(t,x,—D)vw + uvA(t,x,D)w dx 
The main bound for Eq is provided in the next lemma. 

Lemma 5.3. Let 1 < /j, < A. Assume that £g is a Lipschitz function of x with 
|£e — 0\ <C 1 and that a € C 1 5'^ om - Then the trilinear form Eq satisfies the fixed 
time estimate: 

\E (u,v,w)\ < \\u\\ L p 1 \\v\\ L < n \\w\\ L r 1 

+ A _1 ||(^6/ A D)u\\ L p 2 ||u||l«2 I (£e A D)w\\l^ 2 
+ A" 1 !|w||lp2 ||(fr A D)v\\ Lq2 ||(& A D)w\\ L r 2 
+ »\- 2 \mAD)u\\ L v 3 \\(ligAD)v\\ L * 3 \\w\\ L r 3 



(90) 
for all indices 



1 1 1 

1 1 = 1, 1 < Pi, qi,n < oo 

Pt q% n 



and for all functions u, v localized at frequency X in a small angular neighbourhood 
of 8, respectively —6 and all w localized at frequency ll. 

While any choice of LP norms is allowed in the lemma, in order to conclude the 
proof of the estimate for E it suffices to use the set of indices (2, 2, oo). We apply 
the lemma with u = ug 1 v = vg and w — Lwg as above. This yields 



/ 



E (ug,Vg,w)dt 



< IMIl-lHMU^IIwIIl^l-. 



(91) + A^IK^e A D)ug\\ L ^ L 2\\vg\\ Lx , L 2\\(£ e A D)w\\ L 2 Lx , 

+ A^IMUoo^lKfr A D)v e \\ L oo L 2\\(te A D)w\\v> L ~ 
+ li\- 2 \\{ie f\D)ug\\ L ^ L 2\\(ie f\D)v 6 \\ L oov\\w\\v L oa 

Due to the angular localization, the operator (£g AD) yields a factor of ^T^X when 
applied to ug or vg, respectively a factor of a/i when applied to w. Hence we obtain 



/ 



E (ug,Vg,w)dt 
which is acceptable since a' 2 [i> 1. 



IX 2 Ql 2 i i 

< 5 (1 + OlH* + Oi\l2 + l)\\ug\\ x \, a ,e\\vg\\ x x, a ,e\\wg\ 

a z /j, + - 
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Proof of Lemma \5.Sk Since the symbol a is smooth and homogeneous of order 1 
with respect to £, we can use its representation in terms of the spherical harmonics 
and reduce the problem to the case when a has the form 

a(x,Z)=b(x)c{Q 

where b is Lipschitz continuous. 

We denote by £, respectively r\ the frequencies for the ug, respectively vg factors 
in Eq. Then £ and r\ have size A and are in a small angular neighbourhood of 6. 
We expand c around the line generated by £g into a linear term and a quadratic 
error, 

c(0=£(Vc)(&) + £B(£,&)£ 
where B is homogeneous of order —1 with respect to £ and can be chosen so that 

&£(o&) = o, £(&&)& = o 

To see that this is possible we observe that after a rigid rotation we can assume that 
£e = e\. For £ = (1,0) with |0| <C 1 we write the first order Taylor polynomial 
with integral remainder 

c(l, 0) = c(l, 0) + £'c e (1, 0) + ?B(1, 0)0 

= c?1 (i, o) + 0c e (1,0) + 05(i, 0)0 

where -B is given by 

5(1,0)= / (l-/i)Vf,a(l,/i0)* 
Jo 

This extends by homogeneity to all £ in a small angle around 9. 
We represent B as a rapidly convergent sum of terms of the form 

where g is a scalar function which is bounded and smooth on the A scale and F is 
a matrix inheriting the above property of B, 

(92) teFfa) = 0, F(&)& = 

So we have 

c(0 = £(Vc)(&) + A- 1 E ^(&)&(0 
Then we obtain the rapidly convergent series representation 

c(0-cfa) = ^-^(V^^^ + A-^^-^^O 
+ A- 1 ^^)^-^^) 

+ a- 2 e^(&)£(£-^(0%) 

where /i and fc are smooth and bounded on the A dyadic scale. 

We use this representation for the first two components in Eq. The contribution 
of the first term above cancels the principal part of the third component in Eq. 
We retain the other three terms though, therefore this yields the following rapidly 
convergent series representation for Eq: 

E (u, v,w)= ( uvwD(b(Vc)(&))dx + ^ Eq + E E o + E E o 
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The first term is easily estimated since 6(Vc)(£g) is Lipschitz continuous. The first 
summand has the form 

E\ = A" 1 J F{ig)D{Dg{D)uv) wdx 

= -A" 1 J DF{£ e )Dg(D)v,vw + F(£ g )Dg(D)uvDw dx 

In the first term F{£g) is Lipschitz in x and the u derivative yields a factor of A. 
For the second term on the other hand we use (1521 to estimate 



\Dg(D)uF£ e )Dw\ < \fa A D)g(D)u\\^g A D)w\ 

Commuting g(D) with (£g A D) we get 

\{ZeAD)g{D)u\ < \g(D)(t; g A D)u\ + \[g(D),Z e A D]u\ 

with the commutator [g(D),£g A D] bounded in all L p spaces. Hence 

l^ol £ IMU" IMU« \\w\\li + A -1 IK^e A D)u\\ LP2 \\v\\ L , 2 ||(& A D)w\\ L n 

The second summand of Eq is similar but with the roles of u and v reversed. 
Finally, 

Eq = \- 2 J Ffo) D(Dh{D)uDk(D)v) wdx 

= -A" 2 J DF{£_ ) Dh(D)u Dk(D)v w + F(£g) Dh{D)u Dk(D)v Dw dx 

where the matrix F(£$) is paired with the u and v derivatives. In the first term the 
two derivatives on u and v yield a A 2 factor. In the second term we use as before 
([92]) and commute out the h{D) and k(D) multipliers. We obtain 

l^ol £ \\u\\lpi\\v\\lii\\w\\l-i + fJ.\- 2 \\{& A D)u\\ LP3 \\{£ 9 A D)v\\ Lq3 \\w\\ L r 3 

Summing up the results we get the conclusion of the Lemma. 

□ 

Case 1, term III: This has the form 

111= I ^2 ^2 ^i,'x(^ x ' D ) S ^ u 4>Z s a x ( t ,x, D )S\ v 4>t,'"^ X {t,x,D)S^wdxdt 
J i eeOc 

In this case the summation with respect to 9 is accomplished by (I77p , while for the 
a summation we simply accept a In /j, loss. Fixing a and 9 we set 

u e = 4>t,\(t,x, D )S\u, Vg = (fTfixfaxjD^xv, w% = 4>j^ x {t,x,D)S fl w. 

and repeat the analysis for Case 1, term II. The angular localization of Ug and Vg 
is not used in the bounds for the first four terms in therefore that part of 

the argument rests unchanged. The same applies to the bound for the fixed time 
integral in ([53]). 

It remains to consider the bound for E(v,g,Vg,w). The a localization angle for 
Wg is now a/i~ 1 A, therefore part (b) of Lemma 15.21 gives 



w\\x+ < -0T2 II ^9 \\x 



fi 
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This is stronger than in the previous case because it gives a high frequency gain. 
Now we are able to use Lemma 15.31 with exponents (3, 2, 6) to obtain 

E (u%,V$,w)dt\ < ||^|| L 2 L 3||^|| L= o L 2||?i|| i 2 L 6 

+ A -1 ||(^ e A £>K|U 2 l3|KIIl^||(& A D)w\\ L 2 L e 

+ A- 1 ||ug , || L » z »||(& A D)v?\\ L ~ L 4fa A D)w\\ L 2 L e 

+ n\- 2 \m A D)u$\\v L s\\fa A £>K|| L oo i2 |H| i2i a 

Due to the angular localization on the a scale for Ug and Vq , respectively on the 
a/i _1 A scale for Wg, all (£g A D) operators above yield aA factors. Hence, taking 
advantage of the Strichartz estimates, we obtain 

E (u%,v%,w)dt\ < ~^a 2 X \&^u$\\ x *. a ,4i%\\ x *, a ,.\\u%\\ „^, e 
a A X z + - f ' 

= A-A/xv|| u «|| x x, a ,,\\vf\\ x x,*,4v%\\ 

A + 

which is satisfactory since A > /x. 

We conclude this case with two remarks. First, in this context the proof of 
Lemma [5.3l is somewhat of an overkill. In fact, it would suffice to linearize separately 
a{t, x,£) and a(t, x, rj) around £g and use the fact that the symbol a(t, x, £) — a(t, x, £) 
has size a 2 X at frequency A in H a S a {9). Secondly, the endpoint Strichartz estimate 
is only used here for convenience; there is some flexibility in choosing the indices. 

Case 1, term IV. This has the form 

IV = I ^2 ^ <t>e'\( t > x > D ) s * u <t>-'e 1 .\(t> x -> D )Sx v < t > 6.'u M X {t 7 x,D)S^wdxdt 
■> i eeo a 

Again the summation with respect to 9 is accomplished by ([77]) . while for the a 
summation we simply accept a In /i loss. This term is better behaved because the 
symbol 

vanishes on H = {£ + r\ + £ = 0}. Precisely, in the support of the above symbol we 
have 

lf|«A, \tl\$\<a\, H«A, \ V A$\ttCa\, |C| « A, |C A $\ < aX. 
This leads to 

(93) |(£ + »7 + 0A#|»CaA 

This can be taken advantage of in a direct computation in the above formula. 
Including the dyadic frequency localizations into the <f>'s, each term in TV has the 
integral representation 

x, OHO fcftfe x, T))m 4r~ lA (*' C ^( C ) e^^+O d^drjdCdxdt 
Defining the spatial elliptic operator F with symbol 

we have 

F(t,x,D x )e lx ^ +T ' +c} = |(C + 7? + C) A$\ 2N e lx ^ +7 i+V 
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Hence integration by parts in the above formula leads to 

ip(t, x, f , T), C)w(0 v(v) w(()e lx< t+v+0 d£drjd(dxdt 
where the new symbol tp is 



iP(t,x,£,ri,0 = F*(t,x,D x ) 



|(£ + »7 + C)A#l 



a\2N 



In the support of the numerator the bound (|93| holds. Hence separating the vari- 
ables we can represent the denominator as a rapidly convergent series with terms 

(a\r 2N X <ax(H A $) X Cax(v A &)x< a x(C A €$) 

where each of the x's above is a unit bump function on the aA scale. Thus they 
can be included in the corresponding <fi factors. Due to the S(g a ) regularity of the 
4> factors, each derivative ^°Afl applied to them yields an a -1 factor. Thus ip is 
represented as a rapidly convergent series of products of the form 

where the ip factors have the same support and regularity as the corresponding </>'s. 
The integral above is similarly represented as a rapidly convergent series with terms 
of the form 

{a 2 X)- 2N J tpgf(t, x, D)S x u C£ A (t, x, D)S x v ^^{t, x, D)S^w dxdt 

Since a > fjT ' , the factor in front of the above integral allows us to exchange low 
frequencies for high frequencies. This suffices in order to bound the last integral 
using Strichartz estimates. 
Case 1, term V 

This is similar to Case 1, term IV. This time in the support of the symbol 

<t>X t (*> x > o (*> x ' <t>i;» (*« x > 

we have 

Id ~ A, |e A ^1 < «A, |r,|«A, \r)A$\*Ca\, \(\ « A, |C A $\ « C(3\. 
Hence 

|(£ + »7 + C)A#|«CaA 
therefore the symbol above is supported at distance (3X from the diagonal H . Hence 
integrating by parts as in the previous case we gain arbitrary powers of (a/3A) _1 . 
Then we can close the argument using Strichartz type estimates. 

Case 2, 1 < d < fi. This requires only minor changes, which we describe in 
what follows. We still consider the five terms in the trilinear decomposition (|5Tj) . 
but we replace the smallest localization angle yT^ by d^pT^ . 

Case 2, term I. Here we need (|77|) to sum expressions of the form 

A+ A d2 '' 2 (t, x,D)S x u (f>Z'e,\ >i 2 (t,x, D )S\ v < l ) e^ 2>1 2 (t,x,D)S l _ l wdx 

over € O^i _ i . Each term is bounded by combining the energy estimate for the 
first factor, the L A L 2 bound for the second and (|73| for the third. 
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Case 2, term II. Here we use (|77|) for the summation of expressions of the form 

r l _ l l _ l 



J i _ i 

over 9 S O 1 1 . We use the same operator L, the same function w and the same 

trilincar form E. In (|88[) the first, second and fourth terms are estimated in the 
same way, but using the L A L 2 bound for the second factor. In the third term we 
lose a power of d, 



uo(D t — A(t, x, —D))vgwdx 



< 



< 



< 



\ue\\L°°L'\\(Dt ~ A(t,x,-D))v g \\ L 2\\w\\ 
\ug\\ x+ di\\vg\\x_^ d -^-aifi^\\w^\\ x ^, 

3 



a 2 fi 



But this is still acceptable due to the reduced range for a, namely a 2 [i > d. 

In the expression ([55)1 there is a d* loss in the L 2 bound for vg, but this is 
compensated for by the previously unused (a 2 /x) _ * factor in the pointwise bound 
for w. 

Finally, for the Eq bounds we reuse ([91]) but with all the vg factors estimated in 
L 2 . This produces an extra d~4 gain. On the other hand, the angular localization 
for ug and vg is worse. Precisely, the operator (£g A D) yields a factor of d? fi~iA 
when applied to ug or vg, respectively a factor of a/i when applied to w. Hence we 
obtain 

2 I 

U2(J2 1111 

Eo(ug, vg,w)dt < — (1+d? an? +d? an? +d)\\ug\\ x \ ,e,c.\\vg\\ x x.e,c.\\wg\\ x ^,e,, 

dia 2 /j, + _ + 

This is still acceptable since a 2 /i > d. 

Case 2, term III. Compared to the similar argument in Case 1, the following 
modifications are required: 

(i) The third term in ([55]) is treated as in Case 2, term II. 

(ii) In the bound for Eo, the L°°L 2 norms are replaced by L A L 2 in all the vg 
factors. 

Case 2, terms IV, V. These are identical to Case 1. 

□ 
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